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It is shown that the hadron production in high energy pp and pp collisions, 
calculated by assuming that particles originate in hadron gas fireballs at thermal 
and partial chemical equilibrium, agrees very well with the data. The temperature 
of the hadron gas fireballs, determined by fitting hadron abundances, does not seem 
to depend on the centre of mass energy, having a nearly constant value of about 170 
MeV. This value is in agreement with that obtained in e + e~ collisions and supports 
a universal hadronization mechanism in all kinds of reactions consisting in a parton- 
hadron transition at critical values of temperature and pressure. 
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1 Introduction 

The thermodynamic approach to hadron production in hadronic collisions was originally intro- 
duced by Hagedorn |1[ about thirty years ago. The most important phenomenological indication 
of thermal multihadron production in high energy reactions was found in the universal slope 



of the transverse mass (i.e. my = Jpj> + m 2 ) spectra ||, where transverse means orthogonal 
to the beam line. This kind of signature of a hadron gas formation is nowadays extensively 
used in heavy ions reactions, although it has been realized that transverse collective motion of 
the hadron gas may significantly distort the basic thermal m^-spectrum J5|, thus complicating 
the extraction of the temperature. A much better probe of the existence of locally thermalized 
sources in hadronic collisions is the overall production rate of individual hadron species which, 
being a Lorentz-invariant quantity, is not affected by local collective motions of the hadron 
gas. However, the analysis of hadron production rates with the thermodynamical ansatz im- 
plies that inter-species chemical equilibrium is attained, which is a much tighter requirement 
than that of thermal-kinetic intra-species equilibrium assumed in the analysis of ttit spectra Q. 
Chemical equilibrium thus usually implies also thermal kinetic equilibrium. For this reason we 
focus our attention in this paper on the analysis of hadron abundances and the question of 
chemical equilibrium, leaving the analysis of m^-spectra (with its possible complications due 
to collective dynamical effects) to a separate publication. 

The smallness of the collision systems studied here requires appropriate theoretical tools: in 
order to properly compare theoretical predicted multiplicities to experimental ones, the use of 
statistical mechanics in its canonical form is mandatory, that means exact quantum numbers 
conservation is required, unlike in the grand-canonical formalism ||. It will be shown indeed 
that particle average particle multiplicities in small systems are heavily affected by conservation 
laws well beyond what the use of chemical potentials predicts (this was previously observed 
in a similar canonical thermodynamic analysis of pp annihilation at rest ||). However, in the 
high multiplicity (or large volume) limit the grand-canonical formalism recovers its validity. 
This paper generalizes the thermodynamical model introduced in ref. [0] for e + e~ collisions by 
releasing some assumptions which were made there; calculations are performed with a larger 
symmetry group (actually by also taking into account the conservation of the electric charge). 
Moreover, formulae of global correlations between different particles species are provided, and 
a comparison with data is made in this regard as well. 

2 The model 

In refs. 0,|| a thermodynamical model of hadron production in e + e~ collisions was developed 
on the basis of the following assumption: the hadronic jets observed in the final state of a 
e + e~ — > qq event must be identified with hadron gas phases having a collective motion. This 
identification is valid at the decoupling time, when hadrons stop interacting after their forma- 
tion and (possibly) a short expansion {freeze-out). Throughout this paper we will refer to such 
hadron gas phases with a collective motion as fireballs, following refs. fl],|2[]. Since most events in 
a e + e~ — > qq reaction are two-jet events, it was assumed that two fireballs are formed and that 
their internal properties, namely quantum numbers, are related to those of the corresponding 
primary quarks. In the so-called correlated jet scheme correlations between the quantum num- 
bers of the two fireballs were allowed beyond the simple correspondance between the fireball 
and the parent quark quantum numbers. This scheme turned out to be in better agreement 
with the data than a correlation- free scheme (7[. 



The more complicated structure of a hadronic collision does not allow a straightforward exten- 
sion of this model. If the assumption of hadron gas fireballs is maintained, the possibility of 
an arbitrary number of fireballs with an arbitrary configuration of quantum numbers should be 
taken into account |§. To be specific, let us define a vector Q = (Q,N,S,C,B) with integer 
components equal to the electric charge, baryon number, strangeness, charm and beauty respec- 
tively. We assume that the final state of a pp or a pp interaction consists of a set of N fireballs, 
each with its own four-vector /% = Mj/Tj, where Tj is the temperature and Ui = (ji, ficii) is the 
four- velocity flQ| , quantum numbers Q° and volume in the rest frame V^. The quantum vectors 



Q° must fulfill the overall conservation constraint Y^=\ Q° — Q° where Q° is the vector of the 
initial quantum numbers, that is Q° = (2, 2, 0, 0, 0) in a pp collision and Q° = (0, 0, 0, 0, 0) in 
a pp collision. 
The invariant partition function of a single fireball is, by definition: 

S<(Qi°) = E e- ft '% i; , Q o , (1) 

states 

where Pi is its total four- momentum. The factor 8n i qo is the usual Kronecker tensor, which 
forces the sum to be performed only over the fireball states whose quantum numbers Qj are 
equal to the particular set Q°. It is worth emphasizing that this partition function corresponds 
to the canonical ensemble of statistical mechanics since only the states fulfilling a fixed chemical 
requirement, as expressed by the factor 5q.qo, are involved in the sum (1). 
By using the integral representation of ^q^q?: 

r2Trr2Trr2nr2nr2n 



Eq. (1) becomes: 
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*(<$) = EtA^ P ■ ■ A^ e-A'^VWHW* (3) 

states( 27r ) 5 JO J 

This equation could also have been derived from the general expression of partition function 
of systems with internal symmetry ||TT],|12| by requiring a U(l) 5 symmetry group, each U(l) 



corresponding to a conserved quantum number; that was the procedure taken in ref. 0. 
The sum over states in Eq. 3 can be worked out quite straightforwardly for a hadron gas of 
Nb boson species and Np fermion species. A state is specified by a set of occupation numbers 
{rij t k} for each phase space cell k and for each particle species j. Since P, = J2j,kPknj t k and 
Qi = Y*j,k ( \j n j,k, where q^- = (Qj,Nj,Sj,Cj,Bj) is the quantum numbers vector associated to 
the j th particle species, the partition function (3) reads, after summing over states: 

i Nb Np 

ZM) = -— b Jd^ e iQ ^expE£ log(l - e-*™-***)- 1 + £E log(l + e^"^)] . 

\ / j=l k j=l k 

(4) 
The last expression of the partition function is manifestly Lorentz-invariant because the sum 
over phase space is a Lorentz-invariant operation which can be performed in any frame. The 
most suitable one is the fireball rest frame, where the four-vector $ reduces to: 

A = (^,0,0,0) (5) 



Ti being the temperature of the fireball. Moreover, the sum over phase space cells in Eq. (4) 
can be turned into an integration over momentum space going to the continuum limit: 

E-W + Dp^/dV, (6) 

where V is the fireball volume and Jj the spin of the j th hadron. As in previous studies on 
e + e~ collisions 0] and heavy ions collisions JI3], we supplement the ordinary statistical me- 



chanics formalism with a strangeness suppression factor 7 S accounting for a partial strangeness 
phase space saturation^]; actually the Boltzmann factor e _/3 ' Pfc of any hadron species containing 
s strange valence quarks or anti-quarks is multiplied by 7*. With the transformation (6) and 
choosing the fireball rest frame to perform the integration, the sum over phase space in Eq. (4) 
becomes: 



sjg-ft-Pfe-iqj"*^ 1 



£log(l± 7 , 

k 

{2 %y Vi I ^ bg (1 ± 7>-V^/^<^)±i = V . Fj{Ti , 7s , 0) , (7) 

where the upper sign is for fermions, the lower for bosons and Vi is the fireball volume in its rest 
frame; the function Fj(Ti,^y s , <fr) is a shorthand notation of the momentum integral in Eq. (7). 
Hence, the partition function (4) can be written: 

ZM) = j^yj d 5 0e iQ ^exp[V 4 ^F,(T 4 ,7 s ,0)] . (8) 

The mean number (n.,), of the j th particle species in the i th fireball can be derived from Z(Q®) 



by multiplying the Boltzmann factor exp {—Jp 2 + m 2 /T), in the function Fj in Eq. (8) by a 
fictitious fugacity \j and taking the derivative of logZj(Q°, A^) with respect to Xj at A^ = 1: 

(nj)i =^-log^(Q°,A,)| Aj=i . (9) 

The partition function Zj(Qj,Aj) supplemented with the A^ factor is still a Lorentz-invariant 
quantity and so is the mean number (nj)i. From a more physical point of view, this means that 
the average multiplicity of any hadron does not depend on fireball collective motion, unlike its 
mean number in a particular momentum state. 

The overall average multiplicity of the j th hadron, for a set of iV fireballs in a certain quantum 
configuration {Q?, . . . , 0,%} is the sum of all mean numbers of that hadron in each fireball: 

N a a N 











<*i> =^d\~ log Z * (Q * ' Aj) a,=i = ox~ log n z «(Q*> A ^) a,=i • ( 10 ) 

In general, as the quantum number configurations may fluctuate, hadron production should be 
further averaged over all possible fireballs configurations Q°, . . . , Q^ fulfilling the constraint 
Y<iLi Qi — Q°- To this end, suitable weights w(Q?, . . . , Q^r), representing the probability of 
configuration {Q?, . . . , Q%} to occur for a set of N fireballs, must be introduced. Basic features 
of those weights are: 

2 Possible charm and beauty suppression parameters 7 C and 75 are unobservable, see also Appendix C. 



N 

w(q 1 ,...,q°n) = o if £q?^q°, 

I>(Q?, ■ ■ ■ , QSr) = i ■ (11) 

Q?,..,QS, 



For the overall average multiplicity of hadron j we get: 

d 



N 



<K» = E ^(Q?> • • • . Q^)^ log n ZiiQl a,) x =1 . (12) 



Q?,-,QS- 3 i ] 



A,=l 



*N 



There are infinitely many possible choices of the weights w(Q ( j ) , • • • , Q,%), all of them equally 
legitimate. However, one of them is the most pertinent from the statistical mechanics point of 
view, namely: 



*(<& ■■■><&) = f f u :; ov (13) 



^Q°,Q°niIi^(Q°) 

EQO,...,Qo5 Si QO,QO Hill ^i(Qi 

It can be shown indeed that this choice corresponds to the minimal deviation from statistical 
equilibrium of the system as a whole. In fact, putting weights (13) in the Eq. (12), one obtains: 



F) N 

M = or log S **q?,q° II ^(Q. ? . a,o| a _ x . (14) 

3 qO qO j=l J 

This means that the average multiplicity of any hadron can be derived from the following 
function of Q°: 

__ N 

Z(Q°)= E <Wq° I] HQi) . ( 15 ) 

o° O »=1 

with the same recipe given for a single fireball in Eq. (9). By using expression (1) for the 
partition functions Zj(Q°), Eq. (15) becomes: 

N 

m°) = e ^q?,qo n e e- A ^Q0 >Ql . (i6) 

Q°,...,Q° N j=l states; 

Since 

E 5 £;Q?,Q° 5 Q°,Q* = ^Q*,Q° , (17) 

o° o° 

the function (16) can be written as 

Z(Q°) = Y.---Y. *-*'* ■ ■ ■ e- fe ' P ^E iQi , Q o • (18) 

states i statesjM 

This expression demonstrates that Z(Q°) may be properly called the global partition function 
of a system split into N subsystems which are in mutual chemical equilibrium but not in mutual 
thermal and mechanical equilibrium. Indeed it is a Lorentz-invariant quantity and, in case of 
complete equilibrium, i.e. j3i — fa — ■ ■ ■ — (3n = /3, it would reduce to: 



Z{Q°) = E • • • E e-W+'-^fc^.qo = E e-^ Q ,qo , (19) 

states i states^ states 

which is the basic definition of the partition function. 

To summarize, the choice of weights (13) allows the construction of a system which is out of 
equilibrium only by virtue of its subdivision into several parts having different temperatures 
and velocities. Another very important consequence of that choice is the following: if we assume 
that the freeze-out temperature of the various fireballs is constant, that is 7\ = . . . = T/v = T, 
and that the strangeness suppression factor 7 S is constant too, then the global partition function 
(18) has the following expression: 

Z(Q°) = ^ F /dVe 1Q,, -*fflq)[(E t V i )X;F i (r,7.^)]. (20) 

Here the V^'s are the fireball volumes in their own rest frames; a proof of (20) [0] is given in 
Appendix A. Eq. (20) demonstrates that the global partition function has the same functional 
form (3), (4), (8) as the partition function of a single fireball, once the volume Vi is replaced 
by the global volume V = J2f=i Vi- Note that the global volume absorbs any dependence of the 
global partition function (20) on the number of fireballs N. Thus, possible variations of the 
number N and the size Vi of fireballs on an event by event basis can be turned into fluctuations 
of the global volume. In the remainder of this Section and in Sects. 3, 4 we will ignore these 
fluctuations; in Sect. 5 it will be shown that they do not affect any of the following results on 
the average hadron multiplicities. 

The average multiplicity of the j th hadron can be determined with the formulae (14)-(15), by 
using expression (20) for the function Z(Q°): 



((nj)) = 7^5 / ^ e^exp^E^T* 

(21) 



(27T) 

(2Jj + l)V r d 3 p 



(2tt) 3 J 7s S J exp Up 2 + m?/T + ity • 0) ± 1 



where the upper sign is for fermions and the lower for bosons. This formula can be written in 
a more compact form as a series: 

M = Efrirsr^w z(Q °~o? j) > ( 22 ) 



n=\ 



Z(Q°) 



where the functions Zj^ are defined as: 



V r i VT nm ■ 

z j{n) = (2 J, + 1) ^ J d 3 p exp (-ny/p* + mj/T) = (2 J, + 1) Z-. m ) K 2 (-^) . (23) 

K 2 is the McDonald function of order 2. Eq. (22) is the final expression for the average 
multiplicity of hadrons at freeze-out. Accordingly, the production rate of a hadron species 
depends only on its spin, mass, quantum numbers and strange quark content. 
The chemical factors Z(Q° — ncy)/Z(Q°) in Eq. (22) are a typical feature of the canonical 
approach due to the requirement of exact conservation of the initial set of quantum numbers. 
These factors suppress or enhance production of particles according to the vicinity of their 



quantum numbers to the initial Q° vector. The behaviour of Z(Q) as a function of electric 
charge, baryon number and strangeness for suitable T, V and 7 S values is shown in Fig. 1; for 
instance, it is evident that the baryon chemical factors Z(0, N, 0, 0, 0)/Z(0, 0, 0, 0, 0) connected 
with an initially neutral system play a major role in determining the baryon multiplicities. The 
ultimate physical reason of "charged" particle (q^ ^ 0) suppression with respect to "neutral" 
ones (qj = 0), in a completely neutral system (Q° = 0), is the necessity, once a "charged" 
particle is created, of a simultaneous creation of an anti-charged particle in order to fulfill the 
conservation laws. In a finite system this pair creation mechanism is the more unlikely the more 
massive is the lightest particle needed to compensate the first particle's quantum numbers. For 
instance, once a baryon is created, at least one anti-nucleon must be generated, which is rather 
unlikely since its mass is much greater than the temperature and the total energy is finite. On 
the other hand, if a non-strange charged meson is generated, just a pion is needed to balance 
the total electric charge; its creation is clearly a less unlikely event with respect to the creation 
of a baryon as the energy to be spent is lower. This argument illustrates why the dependence 
of Z(Q) on the electric charge is much milder that on baryon number and strangeness (see 
Fig. 1). In view of that, the dependence of Z(Q) on electric charge was neglected in the 
previous study on hadron production in e + e~ collisions J7j. These chemical suppression effects 
are not accountable in a grand-canonical framework; in fact, in a completely neutral system, all 
chemical potentials should be set to zero and consequently "charged" particles do not undergo 
any suppression with respect to "neutral" ones. 

A compact analytic expression for the function Z(Q) does not exist. However, an approximation 
of Z(Q) valid for large global volumes (see Appendix B) exists in which chemical factors reduce 
to a product of a chemical-potential-like factor and an additional multivariate gaussian factor 
having no correspondence in the grand-canonical framework. The gaussian factor tends to 1 for 
V —>■ oo proving the equivalence between canonical and grand-canonical approaches for large 
systems. 

The global partition function (18) has to be further modified in pp collisions owing to a major 
effect in such reactions, the leading baryon effect ||14||. Indeed, the sum (18) includes states with 



vanishing net absolute value of baryon number, whereas in pp collisions at least one baryon- 
antibaryon pair is always observed. Hence, the simplest way to account for the leading baryon 
effect is to exclude those states from the sum. Thus, if \N\ = X)i \Ni\ denotes the absolute value 
of the baryon number of the system, the global partition function (18) should be turned into: 

*=£•••£ «-*•* . . . e-^fc.d.oo -£•••£ ^ Pl ■ ■ ■ e"^(Wl»<W ■ (24) 

statesi states^ statesi states^ 

The first term, that we define as Zx(Q ), is equal to the function Z(Q°) in Eqs. (18), (20), 
while the second term is the sum over all states having vanishing net absolute value of baryon 
number. The absolute value of baryon number can be treated as a new independent quantum 
number so that the processing of the partition function described in Eqs. (l)-(3) can be repeated 
for the second term in Eq. (24) with a U(l) 6 symmetry group. Accordingly, this term can be 
naturally denoted by Z 2 (Q°,0), so that Eq. (24) reads: 

Z = Z 1 (Q°)-Z 2 (Q°,0). (25) 

By using the integral representation of 5|at|,o 

W=^jf%M iW (26) 



in the second term of Eq. (24), one gets: 

^ 2 (Q°,0) = ^jd^e^exp^F^T,^)] 

x j<ty ex P E ^|^ /d^log(l + 7ye-^^/ y -^^-^)] , (27) 

where the first sum over j runs over all mesons and the second over all baryons. The average 
multiplicity of any hadron species can be derived from the global partition function (25) with 
the usual prescription: 

(M = £-togZ(X j )^_. (28) 



d\j 



A,=l 



3 Fit procedure and data set 

The model described so far has three free parameters: the temperature T, the global volume 
V and the strangeness suppression parameter j s . They will be determined by a fit to the 
available data on hadron inclusive production at each centre of mass energy. Eq. (22) yields 
the mean number of hadrons emerging directly from the thermal source at freeze-out, the so- 
called primary hadrons [[7|,|I5]|, as a function of the three free parameters. After freeze-out, 
primary hadrons trigger a decay chain process which must be properly taken into account in a 
comparison between model predictions and experimental data, as the latter generally embodies 
both primary hadrons and hadrons generated by heavier particles decays. Therefore, in order 
to calculate overall average multiplicities to be compared with experimental data, the primary 
yield of each hadron species, determined according to Eq. (22) (or (28) for pp collisions) is 
added to the contribution stemming from the decay of heavier hadrons, which is calculated by 



using experimentally known decay modes and branching ratios |16|. [L7| . 

The calculation of the average multiplicity of primaries according to Eq. (22) involves several 
rather complicated five-dimensional integrals which have been calculated numerically after some 
useful approximations, described in the following. Since the temperature is expected to be below 
200 MeV, the primary production rate of all hadrons, except pions, is very well approximated 
by the first term of the series (22): 

«^ £ w 1 - (29) 

where we have put Zj = Zj(i). This approximation corresponds to the Boltzmann limit of Fermi 
and Bose statistics. Actually, for a temperature of 170 MeV, the primary production rate of 
K + , the lightest hadron after pions, differs at most (i.e. without the strangeness suppression 
parameter and the chemical factors which further reduce the contribution of neglected terms) 
by 1.5% from that calculated with Eq. (29), well within usual experimental uncertainties. 
Corresponding Boltzmannian approximations can be made in the function Z(Q), namely 



log (1 ± e - V / P 2 + m "/ T - i ^>)±l ^ g-yV+mf/T-iq^ ^ ^ 

which turns Eq. (20) (for a generic Q) into: 
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Z(Q) ^J^yJ^ eiQ ^ ex P E zjrfe- 1 *-* + J2^y 3 Jd 3 p log (1 - e -y/?*%/T-W)-i] , 

(31) 
where the first sum runs over all hadrons except pions and the second over pions. 
As a further consequence of the expected temperature value, the z functions of all charmed 
and bottomed hadrons are very small: with T = 170 MeV and a primary production rate of K 
mesons of the order of one, as the data states, the z function of the lightest charmed hadron, 
D°, turns out to be ~ 10~ 4 ; chemical factors produce a further suppression of a factor « 10~ 4 . 
Therefore, thermal production of heavy flavoured hadrons can be neglected, as well as their z 
functions in the exponentiated sum in Eq. (31), so that the integration over the variables </>4 
and 4>5 can be performed: 

Z(Q,C,B) ~ * Id^e^expE^e- 1 ^ 

- E 7^3 / d 'P lQ g (1 - e-V^/ r - u u*)] 5 c ,oS B ,o = C(Q) *c,«Aj,o - (32) 



Q and q^ are now three-dimensional vectors consisting of electric charge, baryon number, and 
strangeness; the five-dimensional integrals have been reduced to three-dimensional ones. 
Apart from the hadronization contribution, which is expected to be negligible in this model, 
production of heavy flavoured hadrons in hadronic collisions mainly proceeds from hard per- 
turbative QCD processes of cc and bb pairs creation. The fact that promptly generated heavy 
quarks do not reannihilate into light quarks indicates a strong deviation from statistical equi- 
librium of charm and beauty, much stronger than the strangeness suppression linked with 7 a . 
Nevertheless, it has been found in e + e~ collisions (7] that the relative abundances of charmed 
and bottomed hadrons are in agreement with those predicted by the statistical equilibrium 
assumption, confirming its full validity for light quarks and quantum numbers associated to 
them. The additional source of heavy flavoured hadrons arising from perturbative processes 
can be accounted for by modifying the partition function (31). In particular, the presence of 
one heavy flavoured hadron and one anti-flavoured hadron should be demanded in a fraction 
of events / = cr(pp(p) — > cc)/cr(pp(p)) (or / = cr(pp(p) — > bb)/a(pp(p))) where cr(pp(p)) is 
meant to be the total inelastic or non-single-diffractive cross section. Accordingly, the partition 
function to be used in events with a perturbative cc pair, is, by analogy with Eq. (24)- (25) and 
the leading baryon effect: 

Z =J2---J2 e-^...e-^^ SiQi , Q o 

statesi statesjM 

E • • • E e_/3rPl • • • e-^- Pjv <W iQ o5| C | i0 = Si(Q°) - Z 2 (Q°, 0) , (33) 

states i states^ 

where \C\ is the absolute value of charm. The primary yield of charmed hadrons, calculated 
according to Eq. (28) and partition function (33), is derived in Appendix C. 
A significant production rate of heavy flavoured hadrons might affect light hadrons abundances 
through decay feed-down, so it is important to know how large the fraction / is. Available 



data on charm cross-sections [18| indicate a fraction / « 10 2 -j- 10 3 at centre of mass energies 

9 



Table 1: Values of fitted parameters in pp and pp collisions. The normalization parameter VT 3 is better suited 
than V in the fit because is less correlated to the temperature. The additional errors within brackets have been 
estimated by excluding some data points and repeating the fit. Also quoted is the correlation parameter between 
T and VT 3 . 



V~s (GeV) 


T (MeV) 


VT 3 


7s 




x 2 / dof 


p(T,VT 3 ) 


pp collisions 


19.4-=- 19.6 

23.8 

26.0 

27.4 ~ 27.6 


190.8 ±27.4 
194.4 ± 17.3 
159.0 ±9.5 
169.0 ±2.1 (±3.4) 


5.79 ±3.05 
6.34 ±2.49 
13.36 ± 2.66 
11.04 ±0.69 (±1.4) 


0.463 ± 0.037 
0.460 ± 0.067 
0.570 ± 0.030 
0.510 ±0.011 


(±0.025) 


6.38/4 
2.43/2 
1.86/2 
136.4/27 


-0.999 
-0.936 
-0.993 
-0.972 


pp collisions 


200 
546 

900 


175.4 ± 14.8 
181.7 ± 17.7 
170.2 ± 11.8 


24.26 ± 7.89 
28.5 ± 10.4 
43.2 ± 11.8 


0.537 ±0.066 
0.557 ±0.052 
0.578 ± 0.063 




0.698/2 
3.80/1 

1.79/2 


-0.989 
-0.993 
-0.982 



< 30 GeV and, consequently, much lower values for bottom quark production. Therefore, 
the perturbative production of heavy quarks can be neglected as long as one deals with light 
flavoured hadron production at yfs < 30 GeV. We assume that it may be neglected at any 
centre of mass energy; this point will be discussed in more detail in the next section. 
All light flavoured hadrons and resonances with a mass < 1.7 GeV have been included among 
the primary generated hadron species; the effect of this cut-off on obtained results will be 
discussed in the next section. The mass of resonances with T > 1 MeV has been distributed 
according to a relativistic Breit-Wigner function within ±21" from the central value. The 
7 S strangeness suppression factor has also been applied to neutral mesons such as 0, u, etc. 
according to the their strange valence quark content; mixing angles quoted in ref. [|T(J have 
been used. Once the average multiplicities of the primary hadrons have been calculated as a 
function of the three parameters T, V and 7 S , the decay chain is performed until n, /i, K 1 * 1 , 
K°, A, 2, S 1 * 1 , Q~ or stable particles are reached, in order to match the average multiplicity 
definition in pp and pp collisions experiments. It is worth mentioning that, unlike pp and pp , 
all e + e~colliders experiments also include the decay products of K°, A, H, E^ and Q~ in their 
multiplicity definition. 
Finally, the overall yield is compared with experimental measurements, and the \ 2 : 

X 2 = yXtheOj — expej) 2 /error^ (34) 

i 

is minimized. 

As far as the data set is concerned, we used all available measurements of hadron multiplicities 
in non-single-diffractive pp and inelastic pp collisions down to a centre of mass energy of about 
19 GeV (see Tables 2 and 3), fulfilling the following quality requirements: 

1. the data is the result of an actual experimental measurement and not a derivation based 
on isospin symmetry arguments; indeed, this model predicts slight violations of isospin 
symmetry due to mass differences; 

2. the multiplicity definition is unambiguous, that means it is clear what decay products 
are included in the quoted numbers; actually, all referenced papers take the multiplicity 
definition previously mentioned; 

3. the data is the result of an extrapolation of a spectrum measured over a large kinematical 
region. 

10 



Some referenced papers about pp collisions quote cross sections instead of average multiplici- 
ties. In some cases (e.g. ref. jOJ) both of them are quoted for some particles, which makes it 
possible to obtain the average multiplicity of particles for which only the cross section is given. 
Otherwise, total inelastic pp cross sections have been extracted from other papers. 
Whenever several measurements at the same centre of mass energy have been available, aver- 



ages have been calculated according to a weighting procedure described in ref. [^OJ prescribing 
rescaling of errors to take into account a posteriori correlations and disagreements of experi- 
mental results. 

Since the decay chain is an essential step of the fitting procedure, calculated theoretical multi- 
plicities are affected by experimental uncertainties on masses, widths and branching ratios of 
all involved hadron species. In order to estimate the effect of these uncertainties on the results 
of the fit, a two-step procedure for the fit itself has been adopted: firstly, the fit has been 
performed with a x 2 including only experimental errors and a set of parameters T , Vq, 7 s o has 
been obtained. Then, the various masses, widths and branching ratios have been varied in turn 



by their errors, as quoted in ref. ||16|| , and new theoretical multiplicities calculated, keeping the 
parameters T ,Vq, 7^0 fixed. The differences between old and new theoretical multiplicity values 
have been considered as additional systematic errors to be added in quadrature to experimental 
errors. Finally, the fit has been repeated with a x 2 including overall errors so as to obtain final 
values for model parameters and for theoretical multiplicities. Among the mass, width and 
branching ratio uncertainties, only those producing significant variations of final hadron yields 
(actually more than 130) have been considered. 

4 Results and checks 

The fitted values of the parameters T, V, 7 S at various centre of mass energy points are quoted 
in Table 1 while the fitted values of average multiplicities are quoted in Table 2, 3 along with 
measured average multiplicities and the estimated primary fraction. The fit quality is very 
good at almost all centre of mass energies as demonstrated by the low values of x 2 's and by the 
Figs. 2, 3, 4, 5, 6. Owing to the relatively large value of x 2 a t \/s = 27.4 GeV, variations of 
fitted parameters larger than fit errors must be expected when repeating the fit excluding data 
points with the largest deviations from the theoretical values. Therefore, the fit at yfs = 27 A 
GeV pp collisions has been repeated excluding in turn (A , p°, (p) and (K~, pions), respectively, 
from the data set; the maximum difference between the new and old fit parameters has been 
considered as an additional systematic error and is quoted in Table 1 within brackets. 
The fitted temperatures are compatible with a constant value at freeze-out independently of 
collision energy and kind of reaction (see Fig. 7). On the other hand, 7 S exhibits a very slow 
rise from 20 to 900 GeV (see Fig. 8); its value of ~ 0.5 over the whole explored centre of 
mass energy range proves that complete strangeness equilibrium is not attained. Moreover, the 
temperature value ~ 170 MeV is in good agreement with that found in e + e~ collisions [p|,[33|1 



and in heavy ions collisions P4"| . On the other hand, the global volume does increase as a 
function of centre of mass energy as it is proportional, for nearly constant T and 7 S , to overall 
multiplicity which indeed increases with energy. Its values range from 6.4 fm 3 at ^/s = 19.4 
GeV pp collisions, at a temperature of 191 MeV, up to 67 fm 3 at -/s = 900 GeV pp collisions 
at a temperature of 170 MeV. However, since volume values are strongly correlated to those 
of temperature in the fit, errors turn out to be quite large and fit convergence is slowed down; 
that is the reason why we actually fitted the product VT 3 instead of V alone. 
Once T, V and 7 S are determined by fitting average multiplicities of some hadron species, their 
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Table 2: Measured average multiplicities compared with fitted theoretical values in pp collisions. The first 
quoted error beside measured values is the experimental error, the number within brackets is the error due to 
uncertainty on masses, widths and branching ratios of the various hadrons. Also quoted are the theoretical 
estimates of the fraction of primaries. 



Particles 



Measurement 



Calculated Primary fraction References 



N /5 = 19.4-M9.7GcV 



Neg. charged 

Charged 

n° 

K2 

P° 

A 

A 



2.85 ± 0.040 (±0.063) 2.798 

7.69 ±0.070 (±0.13) 7.620 

3.34 ±0.24 (±0.11) 3.404 0.173 

0.174 ±0.013 (±0.002) 0.160 0.325 

0.33 ±0.06 (±0.025) 0.448 0.568 

0.0977 ±0.0097 (±0.0056) 0.110 0.243 

0.0136 ± 0.0041 (±0.0007) 0.0135 0.238 



[ill 
[21| 



0.369° 



22 

22 



m\ 

|2|, J21 

j 22 1, 1 21 



y/5 = 23.8 GeV 



7T° 


3.42 ±0.62 (±0.12) 


3.908 


0.166 


\m c 


K 


0.22 ±0.025 (±0.003) 


0.198 


0.319 - 


0.362° 27 


K*++K*- 


0.137 ±0.043 (±0.002) 


0.165 


0.658 - 


0.492 d |2JZ| C 


A 


0.11 ±0.02 (±0.007) 


0.126 


0.238 


[21 


A 


0.021 ± 0.004 (±0.001) 


0.0202 


0.233 


H 



yfs = 26.0 GeV 



Neg. charged 
Charged 

K 

A 

A 



3.53 ±0.05 (±0.094) 3.545 

9.06 ±0.09 (±0.18) 9.087 

0.26 ± 0.01 (±0.005) 0.256 0.507 - 0.559° 

0.12 ±0.02 (±0.009) 0.147 0.295 

0.013 ±0.004 (±0.0007) 0.0120 0.292 



[211 

isi 

[2£| 



v / i= 27.4-^27.6 GeV 



7T U 

K+ 
K" 
K° 

V 

P° 

P + _ 

P 

w 

K*+ 

K*" 

K*° 

K*° 

</> 
f 2 (1270) 

P 

P 
A 
A 

S+ 

£- 

A++ 

A° 

A++ 

A° 

s*+ 

s*- 

A(1520) 



4.10 ±0.11 (±0.15) 4.147 0.293 

3.87 ±0.12 (±0.16) 4.197 0.258 

3.34 ±0.08 (±0.12) 3.269 0.223 

0.331 ± 0.016 (±0.007) 0.302 0.484 

0.224 ±0.011 (±0.004) 0.182 0.380 

0.232 ±0.011 (±0.004) 0.232 0.446- 

0.30 ± 0.02 (±0.054) 0.366 0.453 

0.385 ±0.018 (±0.038) 0.543 0.628 

0.552 ± 0.083 (±0.046) 0.601 0.657 

0.355 ± 0.058 (±0.033) 0.421 0.569 

0.390 ±0.024 (±0.002) 0.443 0.665 

0.132 ±0.016 (±0.002) 0.111 0.742 

0.088 ±0.012 (±0.001) 0.0617 0.628 

0.119 ±0.021 (±0.002) 0.0927 0.679 

0.0903 ± 0.016 (±0.001) 0.0708 0.687 

0.019 ± 0.0018 (±0.) 0.0262 1.00 

0.092 ± 0.012 (±0.002) 0.0684 0.845 

1.20 ±0.097 (±0.022) 1.060 0.337 

0.063 ± 0.002 (±0.001) 0.0610 0.283 

0.125 ±0.008 (±0.008) 0.136 0.276 

0.020 ±0.004 (±0.0008) 0.0147 0.273 

0.048 ±0.015 (±0.004) 0.0423 0.688 

0.0128 ± 0.0061 (±0.0032) 0.0310 0.592 

0.218 ±0.0031 (±0.013) 0.250 0.758 

0.141 ± 0.0098 (±0.0089) 0.212 0.714 

0.013 ± 0.0049 (±0.00049) 0.0111 0.548 

0.0336 ± 0.008 (±0.0006) 0.0165 0.697 

0.020 ±0.0025 (±0.0011) 0.0230 1.00 

0.010 ±0.0018 (±0.0007) 0.0139 1.00 

0.017 ±0.0031 (±0.0005) 0.00996 1.00 



0.495° 



a - The ir° multiplicity is defined in this paper as half the photon multiplicity; therefore, the experimental value 

has been fitted to half the number of photons coming not only from tt° but also from r], u) and S° decays. 

b - Primary fraction of K° and K° respectively 

c - This paper quotes only the cross section. The multiplicity has been obtained by using the total inelastic cross 

section of 32.21 mb at y/s = 23.5 GeV quoted in ref. pli| 

d - Primary fraction of K*+ and K* — respectively 

e - Cross-reference to [pl| 



12 



Table 3: Measured average multiplicities compared with fitted theoretical values in pp collisions. The first 
quoted error beside measured values is the experimental error, the number within brackets is the error due to 
the uncertainty on masses, widths and branching ratios of the various hadrons. Also quoted are the theoretical 
estimates of the fraction of primaries. 



Particle 



Measurement 



Calculated Primary fraction References 



v/I = 200 GeV 












Charged 


21.4 ±0.4 (±0.72) 


21.27 






a 


K» 


0.75 ± 0.09 (±0.009) 


0.783 


0.467 


32 




n 


0.75 ±0.1 (±0.05) 


0.794 


0.291 


32 


b 


A 


0.23 ±0.06 (±0.008) 


0.194 


0.263 






B- 


0.015 ± 0.015 (±0.0002) 


0.0123 


0.579 







^/i =546 GeV 



Charged 


29.4 ±0.3 (±0.96) 


29.25 






a 


K° 


1.12 ±0.08 (±0.012) 


1.139 


0.441 


n 




A 


0.265 ± 0.055 (±0.001) 


0.302 


0.252 






~ — 


0.05 ±0.015 (±0) 


0.0228 


0.567 


m 





^s =900 GeV 












Charged 


35.6 ±0.9 (±1.2) 


35.15 






a 


K° 


1.37 ±0.13 (±0.02) 


1.437 


0.497 


32 
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1.0 ±0.2 (±0.09) 


1.188 


0.301 


32 


b 


A 


0.38 ±0.08 (±0.01) 


0.323 


0.269 






B" 


0.035 ±0.02 (±0) 


0.0258 


0.585 







a - The charged track average multiplicity value quoted in this reference has been increased 

by one as leading particles, assumed to be one charged-neutral nucleon-antinuclcon 

pair per event, were excluded. 

b - The neutron average multiplicity quoted in this reference has been increased 

by 0.5 as leading particles, assumed to be one charged-neutral nuclcon-antinucleon 

pair per event, were excluded. 



values can be used to predict average multiplicities of any other species, at a given centre of 
mass energy. 

Since the dependence of the chemical factors on the global volume V is quite mild in the region 
of interest (see Fig. 2), the hadron density mainly depends on the temperature and 7 S (cf. 
Eqs. (22), (29)). Therefore, constant values of temperature and 7 S imply a nearly constant 
hadron density at freeze-out, which turns out to be « 0.4 -f- 0.5 hadrons/fm 3 , as shown in 
Fig. 10, corresponding to a mean distance between hadrons of approximately ~ 1.6 -j- 1.7 fm. 
Unfortunately, due to its dramatic dependence on the temperature, all density values, except 
that at -y/i = 27.4 GeV, are affected by large errors, and thus a definite claim of a constant 
freeze-out density cannot be made. The same statement is true for the pressure, also shown in 
Fig. 10, whose definition is given in Appendix D. 

The physical significance of the results found so far depends on their stability as a function of the 
various approximations and assumptions which have been introduced. First, the temperature 
and 7s values are low enough to justify the use of the Boltzmann limits (29), (30) for all hadrons 
except pions, as explained in Sect. 3. As far as the effect of a cut-off in the hadronic mass 
spectrum goes, the most relevant test proving that our results so far do not depend on it is 
the stability of the number of primary hadrons against changes of the cut-off mass. The fit 
procedure intrinsically attempts to reproduce fixed experimental multiplicities; if the number 
of primary hadrons does not change significantly by repeating the fit with a slightly lower cut- 
off, the production of heavier hadrons excluded by the cut-off must be negligible, in particular 
with regard to its decay contributions to light hadron yields. In this spirit all fits have been 
repeated moving the mass cut-off value from 1.7 down to 1.3 GeV in steps of 0.1 GeV, checking 
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Table 4: Fit results for pp collisions at y/s = 900 GeV with f c = 0.3 compared with those without perturbative 
charm production. 



Parameter 

Temperature (MeV) 
VT 3 

Is 

X 2 /dof 


fc = 0.3 
163.8 ± 10.9 
46.0 ±12.1 
0.571 ± 0.070 
3.09/2 


/c=0 

170.2 ± 11.8 
43.2 ± 11.8 
0.578 ± 0.063 
1.79/2 



the stability of the amount of primary hadrons as well as of the fit parameters. It is worth 
remarking that the number of hadronic states with a mass between 1.7 and 1.6 GeV is 238 out 
of 535 overall, so that their exclusion is really a severe test for the reliability of the final results. 
Figure 11 shows the model parameters and the primary hadrons in pp collisions at yfs = 900 
GeV; above a cut-off of 1.5 GeV the number of primary hadrons settles at an asymptotically 
stable value, whilst the fitted values for T, V, 7 S do not show any particular dependence on 
the cut-off. Therefore, we conclude that the chosen value of 1.7 GeV ensures that the obtained 
results are meaningful. 
As mentioned in Sect. 3, the perturbative production of heavy quarks has been neglected. This 



is legitimate in low energy pp collisions, where it has been actually measured |T1|, but not 
necessarily in yfs = 0(100) GeV pp collisions, where no measurement exists and one has to 
rely on theoretical estimates. In general, the latter predict very low b quark cross sections, but 
a possibly not negligible c quark production. We used the calculations of ref. |53| according to 
which the fraction / of non-single-diffractive events in which cc pairs are produced (see Sect. 
3) rises as a function centre of mass energy. We repeated the fit for pp collisions at y/s = 900, 
where the fraction / is expected to be the largest, by using the upper estimate of a cross section 
a(pp(p) — » cc) ~ 12 mb, corresponding to / ~ 0.3, in order to maximize the effect of charm 
production. The partition function to be used in such events is that in Eq. (33) with a further 
modification according to Eq. (24) to take into account the leading baryon effect. The model 
parameters fitted with / = 0.3 are quoted in Table 4; their variation with respect to / = is 
within fit errors, implying that extra charm production does not affect them significantly. 

5 Fluctuations and correlations 

In the description of the model and the comparison of its predictions with experimental data, 
we tacitly assumed that the parameters T, V and 7 S do not fluctuate on an event by event 
basis. If freeze-out occurs at a fixed hadronic density in all events, as argued in Sect. 4, then it 
is a reasonable ansatz that T and 7 S do not undergo any fluctuation since the density mainly 
depends on those two variables. However, there could still be volume fluctuations due to event 
by event variations of the number and size of the fireballs from which the primary hadrons 
emerge. 

We will now show that, as far as the average hadron multiplicities are concerned, possible 
fluctuations of V can be reabsorbed in a redefinition of volume provided that they are not too 
large. To this end, let us define p(V) as the probability density of picking a volume between V 
and V + dV in a single event. The primary average multiplicity of the j th hadron is then: 

««,» =Jdv P (v) f>ir +1 7r^ ( n) Z{ % Q ^ ] • (35) 



n=l 
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If the volume V fluctuates over a region where the dependence of chemical factors on it is 
mild (i.e. for large volumes, see Fig. 2), they can be taken out of the integral in Eq. (35) 
and evaluated at the mean volume V. In this case, the integrand depends on the volume only 
through the functions Zj( n ) whose dependence on V is linear (see Eq. (23)) and which can be 
re-expressed as 

z j(n) (V, T, 7 S ) = V i j{n) (T, j s ) . (36) 

Then, from Eq. (35), 

«*i» - E(=Fl) n+1 7re;(n)(T, 7.) Z(( | (Q0 ] qj) / dVp(V) V . (37) 

The integral on the right-hand side is the mean volume V. Thus: 

««*» - E(Tl)" +1 l^^n)(V, T, 7s ) Z(( |° (Q0 ] qj) • (38) 

It turns out that the relative hadron abundances do not depend on the volume fluctuations 
since the mean volume V appearing in the above equation (replacing the volume V in Eq. (23)) 
is the same for all species: all results obtained in Sect. 4 are unaffected. On the other hand, if 
the volume fluctuates over a region where the dependence of chemical factors on it is stronger 
(i.e. in the region of small volumes, see Fig. 2) one can show that the leading term of average 
multiplicities is still given by the Eq. (38) and that further corrections are of the order of D 2 /V , 
where D is the dispersion of the distribution p(V) (see Appendix E). Therefore, if D -C V, 
as it is reasonably to be expected, also in this case the calculation of average multiplicities by 
using a single mean global volume would be a very good approximation. 

We emphasized in the Introduction that the average hadron multiplicities are a very useful 
tool to study hadronization because of their independence from collective dynamical effects. 
More generally, since the number of particles is a Lorentz-invariant quantity, this property is 
shared by the entire multiplicity distribution of any hadron species. However, the shape of the 
multiplicity distribution, unlike its mean value, is affected by volume fluctuations since it is 
actually the superposition, weighted with p(V), of many multiplicity distributions, each of them 
associated with a particular volume V, having different mean values and moments. In a previous 
study |y| it has been shown that the charged particle multiplicity distribution in e + e~ collisions 
at y/s = 91.2 GeV, calculated with a fixed volume, provides a fairly good approximation of the 
experimental data, and that remaining discrepancies between the prediction and the data can 
be explained by assuming a superposition of multiplicity distributions with different volumes. 
This superposition effect (also called shoulder effect) has been further investigated in ref. |36| 



Apart from the mean value, which is the first-order moment, the next lowest order moments of 
multiplicity distribution are related to global correlations between particle pairs. Let us first 
derive them for a fixed volume V: let n) k be the number of hadrons j in the k th phase space cell 



for the i th fireball; then, the overall number of hadron j is J2i k n ) k- According to the partition 



j ' k 

A,k n j,k- 

function (18), the probability of picking a set of occupation numbers {n* >fe }, i.e. a state of the 
system, is 



P({n), k }) = \ exp(- £ A ■ n) m ) 5 Q , Q o . (39) 
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The average number of pairs, whose first member belongs to species j and the second to species 
/, is then: 



(M) =T,(i: n)jd (E <*) n{n), k }) (40) 



states i,k i,k 



if j 7^ /, and 



(( nj{nj 2 1} )) = \ E (E4*) (E4* - i)^(KJ) ( 41 ) 

states i,fc j,fc 

if j = I. In both cases the average number of pairs can be obtained from the partition function 
by multiplying all Boltzmann factors exp (— /% ■ Pk(j)) by fictitious fugacities A/s, one for each 
species, and taking the derivative with respect to Xj, Xi at A = 1: 

((njni - k-,(nj + n,)» = (1 - ^) ^t^^T ,_, ■ ( 42 ) 



Z 9A,-9A; 



A=l 



Since the partition function is Lorentz- invariant and so are the parameters Xj, the average 
number of pairs does not depend, as expected, on the collective fireball dynamics. 
In general, one can show that the average number of n-tuples of K particle species, with rt\ 
particles of species ji, n 2 particles of species j 2 , ■ ■ ., nx particles of species jx, is 

1 : 9 " Z (43) 



ni \...n K \ZdXT ...dXl K ' 

This expression proves that Z(Ai, . . . , A^) is proportional to the generating function of the 
multi-species multiplicity distributions. 

The two-particle global correlation can be defined as the ratio between the actual average 
number of pairs and the one that would have been obtained if their production was independent. 
Thus, if j 7^ I: 

Pl = tfht (44) 

As far as identical particles are concerned, if they were independently produced their multi- 
plicity distribution would be Poissonian and therefore the average number of pairs would be 

<K» 2 A so: 

p .. = M^J (45) 

The calculation of the average number of pairs according to Eq. (42) and the partition function 
(20) yields: 



({n j n l --8 jl {n 2 j +n j ))) 



= (1 - ki) E E(Tl) w+m 7r j+m "(^(n)^MT^(n +w >)) Z(Q ° Z Z J 0) ^^ , (46) 

ra=l n=l \^ I 

where the upper sign is for fermions and the lower for bosons. Whereas the term -2j( n )^(m) is 
present for all particles, the term 8jiZj(n+m) is non-zero only for identical particles; it is a further 
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Table 5: Particle-particle correlations in pp collisions at y/s = 26 GeV |37]] obtained by using the total inelastic 
cross section of 32.80 mb quoted in rcf. |g9[. The errors within brackets next to the theoretical predictions are 
due to finite Monte-Carlo statistics. The experimental values of the correlations p have been estimated by 
dividing the average numbers of pairs by the average multiplicities quoted in ref. [[29| , measured in the same 
experiment. Since the correlation between these two measurements is unknown, the relative experimental error 
on p has been assumed to be the same as on ((711712)). The small effect of quantum statistics on the correlations 
of identical particles at the primary hadron level has been neglected, as explained in the text. 



Particles 

K»A 
K«A 
AA 


(niri2) measured 

0.0530 ±0.0055 
0.0451 ± 0.0052 
0.0079 ± 0.0021 
0.0030 ± 0.0009 


(niri2) calculated 

0.0489(±0.0021) 

0.0472(±0.0018) 

0.0041(±0.0006) 

0.0036(±0.0004) 


p measured 

1.57 ±0.16 
1.45 ±0.17 
2.34 ±0.62 
1.92 ±0.58 


p calculated 

1.49(±0.064) 

1.25(±0.048) 

1.34(±0.18) 

2.05(±0.25) 


p for primaries 

1.44(±0.11) 
1.506 
1.296 
2.820 



contribution to correlated particle production due to quantum statistics, the so called Bose- 
Einstein correlations and Fermi-Dirac anticorrelations. If j 7^ / it turns out that, comparing 
Eq. (46) with Eq. (22), ((rijiii)) 7^ ((ni)) ((ni)) (and there is indeed correlated production), unless 
Z(Q° — nqj —mqi) = Z(Q° — nqj)Z(Q° — mqi)/Z(Q°); this is the case if the function Z(Q) is 
an exponential of Q, which occurs only in the grand-canonical regime (see also Appendix B). 
Thus, in the canonical thermodynamical approach, correlated production of particles belonging 
to different species is definitely an effect of conservation laws in a finite system. As long as 
different species are concerned, owing to the temperature values found in the present analysis 
the contribution to the average number of pairs from terms other than n = 1 and m = 1 in 
series (46) is negligible for all hadrons but pions. Therefore: 

which corresponds to the Boltzmann limit, as discussed in Sect. 3. On the other hand, for all 
identical particles but pions we have: 

..nArij-l).. 1 2 2 Z(Q°-2q,) 

«-^-^» - - 2 lP (*? T z m ) " zm . (48) 

In principle, the term Zj(%) stemming from quantum statistics may not be negligible compared 
to Zj even for high mass hadrons, since the ratio 

*_m 1 K 2 (2mj/r) 
z] 2 Zj K 2 ( mj /T) { ] 

may be of the order of 1 if, due to a very small volume, Zj <C 1 is able to compensate the small 

ratio of McDonald functions. In the present analysis the largest value for the ratio (49) for 

heavy hadrons (i.e. excluding pions) which occurs is 0.096 for K + K + production in pp collisions 

at Vs = 19.4 GeV. 

The global correlation between heavy hadron pairs turns out to be, using Eqs. (44)-(48) and 

(29), 

, r(1 , ^(2) ^(Q°-q J -q^(Q ) ( m 

All calculations performed in this Section refer to primary hadrons, which are not observable 
in actual experiments. Since measured correlations may be affected by the decay chain process, 
the given formulae are not directly comparable with experimental data. Therefore, a complete 
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reconstruction of the production process including both the formation and decay of the primary 
hadrons is necessary in order to test the predictive power of the model in this regard. This can 
be done by a Monte-Carlo procedure: by using the model parameters fitted at each centre of 
mass energy point as described in Sect. 4, a set of numbers {n l ^ k } is generated according to the 
probability (39), and subsequently their decays are performed according to the known decay 



modes and branching ratios as quoted in the Particle Data Book ][L6[ . 

A further problem in the comparison with data is the possibility of volume fluctuations. If the 
volume fluctuates from event to event, the average number of heavy hadron pairs should be 
re-expressed as 

iW ~ \h(n) + *i)» * / dVp(V) i^{ m T S fl z m ) ^^q^' ^ • (51) 

According to what has been stated in the beginning of this section, if we take the chemical 
factors out of the integral and write z = V£(T,j s ), Z( 2 ) — V£(2)(T,j s ), we are left in Eq. (49) 
with both a mean volume (multiplying Sji) and a mean squared volume which is equal to the 
squared mean volume only if the dispersion D of the distribution p(V) vanishes. Taking into 
account volume fluctuations, the correlation pji reads (taking the first term of the series in 
Eq. (38)): 

fV%.* ^(2A ^(Q° - q, - q^(Q°) ( , D 2 z m \ Z(Q - qj - qi )Z(Q ) 

Pjl \V 2 T ]l z) ) Z(Q0 - q^CQO - q,) [ L+ v 2T * z] ) Z(Q° - qj )Z(Q* - q,) ' 

(52) 
The correlation between different particle species then increases by a factor (1 + D 2 /V ) with 
respect to the non-fluctuation case even neglecting the dependence of chemical factors on the 
volume; however, if D <C V this a small effect. 
We compared Monte-Carlo simulated correlations with a set of particle-particle correlations 



measured in pp collisions at y^ = 26.0 GeV [37]; the results are shown in Table 5. The corre 



lations have been predicted by using the model parameters T, V, and 7 S fitted at that centre 
of mass energy quoted in Table 1. We also quote the correlations at the primary hadron level 
which were calculated with Eqs. (47) and (48), taking into account that K° is a mixed particle- 
antiparticle state, and, for the K°K° correlation, with the same Monte-Carlo technique used 
for the final particles. The effect of Bose-Einstein correlations in the global correlated produc- 
tion of K°K°, estimated to be < 1.7% according to formulae (48), (49) and taking mixing into 
account, has been neglected. The comparison between primary and final correlations indicates 
that, in general, they are slightly diluted by the decay chain process. 

Generally, the agreement between the predictions and the data is good. In trusting this ap- 
proach one is led to the conclusion that volume fluctuations are small enough to be hidden in 
the experimental errors. This fact should be confirmed by a more detailed study of charged 
particle multiplicity distributions. 

6 Conclusions 

A detailed analysis of hadron abundances in pp and pp collisions over a large range of centre of 
mass energies (from 20 to 900 GeV) has demonstrated a stunning ability of the thermodynamic 
model to reproduce accurately all available experimental data on hadron production in high 
energy collisions between elementary hadrons. Key elements for the success of this approach are 
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the use of the canonical formalism of statistical mechanics, ensuring the exact implementation 
of quantum number conservation, and the introduction of a supplementary parameter 7 S to 
account for incomplete saturation of strange particle phase space. 

The remarkable agreement of the data with such a purely statistical approach, which uses only 
three free parameters, has important implications. Firstly, it indicates that hadron production 
in elementary high energy collisions is dominated by phase space rather than by microscopic 
dynamics; during hadronization of the prehadronic matter formed in the collision, the hadronic 
phase space is filled according to the law of maximal entropy, with minimal additional (i.e. 
dynamical) information. The only dynamics visible in the final state is the collective motion 
of the hadron gas fireballs which reflects the underlying hard parton kinematics. Secondly, 
the observed universality of the freeze-out temperature independent of the collision energy and 
collision system suggests that hadronization cannot occur before the parameters of prehadronic 
matter, like energy density or pressure, have dropped below critical values corresponding to a 
temperature of around 170 MeV in an (partially) equilibrated hadron gas. Hadronization at 
larger energy densities or pressures is inhibited by the "absence" of a hadronic phase space: 
according to lattice QCD calculations, the most likely (i.e. maximum entropy) state at higher 
energy densities is a colour deconfined quark-gluon plasma in which hadrons do not exist as 
stable degrees of freedom. Therefore, this analysis indicates that the value of critical transition 
temperature is T cr i t ~ 170 MeV. This agrees with the limiting ("Hagedorn") temperature |l],[| 
for an equilibrated hadron gas and with lattice QCD results [J38H . 

The phase-space dominance in the hadronization process can be understood by the non- 
perturbative nature of the strong interaction forces in this energy density domain: within each 
fireball, many different processes and channels contribute to the formation of soft hadrons, 
resulting locally in equal transition probabilities for all hadronic states in phase space. The 
value of temperature reflects the hadronic energy density or pressure at its critical value where 
hadron production occurs while the only other parameter entering the observed hadron spectra 
is the collective motion relative to the observer. 

The only deviation from this picture of complete phase-space dominance in hadronization re- 
sides in the incomplete saturation of strange particle phase space, i.e. 7 S ~ 0.5: the final 
hadronic state seems to "remember" that there were no strange quarks in the initial state, and, 
in spite of their non-perturbative nature and the many possible dynamical channels, strong 
interactions in the pre-hadronic stage do not manage to wipe out completely the asymmetry 
between strange quark and light quark abundances. Strangeness suppression, as well as the 
survival of perturbatively created cc and bb pairs, are thus the only trace to strong interaction 
dynamics before hadronization. The systematics of the observations suggest that these non- 
equilibrium effects are mainly related to quark mass thresholds. Similar analyses of hadron 
abundances in nuclear collisions suggest that the strangeness suppression disappears in larger 



collision systems with larger lifetimes prior to hadron freeze-out [p9 |. Note that also here, in 
hadronic collisions, the slight increase of 7 S with centre of mass energy in hadronic collisions is 
connected with a systematic increase of the fitted fireball volumes at freeze-out (see Table 1); 
this goes in the same direction. The increase of freeze-out volume with rising centre of mass 
energy may imply a corresponding increase of the initial (prehadronic) energy density since the 
final energy density of the hadronic state is limited by the observed constant temperature. 
It has been shown that not only the average single hadron abundances, but also the particle- 
particle correlations measured in pp collisions agree well with the thermal predictions. It should 
be stressed that the requirement of exact quantum number conservation yields major effects on 
both the average hadron multiplicities and the correlations, and that for elementary high en- 
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ergy collisions a thermal description in the grand-canonical framework would not have worked 
so well. As noted by Hagedorn many years ago ||, J5| and extensively discussed in Sect. 2, 
small fireballs volume result in the suppression of strange relative to non-strange hadrons even 
for 7 S = 1 (i.e. without strange phase space suppression), when one compares the canonical 
with the grand-canonical approach, due to the need to create strange particles always in pairs. 
A strangeness suppression 7 S ~ 0.5 in the canonical approach, as extracted here from the pp 
and pp data, thus corresponds to a seemingly much stronger suppression of 7 S ~ 0.2 within a 
grand-canonical approach J5JJ. Vice- versa, it should be emphasized that even a constant value 
of 7 S may imply a strong relative enhancement of strange particles going from the small volumes 
of elementary hadron collisions to possible large volumes in nuclear collisions. The frequently 
discussed "strangeness enhancement" in nuclear collisions thus really consists of two compo- 
nents: (i) the removal of the suppression (at constant j s ) arising from the need to conserve 
exactly strangeness in a small collision volume, and on top of that (ii) additionally a possibly 
larger value of 7 S 0,[39|. Both of these effects are dynamically non-trivial. 
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7 Appendix 

A Proof of Equation (20) 

We want to prove that the global partition function (18) can be expressed by Eq. (20) if the 
temperatures and the strangeness suppression factors 7 S of the various fireballs are constant. 
Let rij k be the number of the j th hadron species in the k th phase space cell of the i th fireball. 
Then: 

E Qi = E n ),k^ 3 

i i,j,k 

Pi = J2<kPk- (53) 

3,k 

By using Eq. (2) and putting Eq. (53) into Eq. (18) the following expression of global partition 
function is obtained: 

1 r N 

z (Q ) = 7o^5 / d 5 0e iQO -*n E exp [-]>>• n> fc - in)^ ■ <j>] . (54) 

V / i=\ statesi j,fc 

After summing over states and inserting the strangeness suppression factor 7 S , the Eq. (54) 
becomes: 

1 r N 

Z(Q°) = — - / d 5 e ! ^n«P EE log(l± T ^e-ft^- i ^^ ±1 ] , (55) 

y A7r ) J i=i j k 

where the upper sign is for fermions and the lower for bosons. 

Once the transformation (6) has been applied in Eq. (55), one is left with phase space integrals 

that may be performed in the rest frame of each fireball, in the very same way as in Eq. (7): 

Z ^ = TVTs I d ^ e i( ^nexp^E^( T ^,0)] • (56) 

\ A7r ) J i=i j 

lfT 1 = ... = T N = T and 7sl = . . . = 7sAr = 7s , then: 

Z(Q°) = j^y 5 J d 5 e^^expf^^^F,^,^,^] . (57) 

which is precisely the Eq. (20). 

B Approximation of the function Z(Q) for large systems 

We look for an approximated expression of the function Z(Q) for large values of particle mul- 
tiplicity, namely for large values of volume V. In the following calculations heavy flavoured 
particles, whose multiplicities are orders of magnitude below light flavoured ones, at temper- 
atures T = 0(100) MeV, are completely neglected. This means that we are dealing with a 
function Z(Q) as in Eq. (32) and that vectors Q and q^ are henceforth meant to be three- 
dimensional with components electric charge, baryon number and strangeness respectively. 
Let us define: 
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M = E 2j ^ j d 3 P log (1 ± tfe-y/F^/r-***)* , (58) 

where the upper sign is for fermions and the lower is for bosons. By using this definition, the 
function C(Q) i n Eq. (32) can be written: 

C(Q) = j^y 3 J d 3 e i( ^exp [Vf{<t>)\ . (59) 

The logarithm in the function /(</>) in Eq. (58) can be expanded in a series: 

oo z-j-i \n+l 

log (1 ± 7 s ^eV p2+m ? /T e~ iq ^' </> ) ±1 = V K ' 7 ns ie" n v p2+m ? /T e~ ni ^^ (60) 

n=l n 

The sum labelled by index j in the function /(0) obviously runs over all particles and anti- 
particles species. However, in order to develop calculations, it is advantegeous to group particles 
and corresponding anti-particles terms together in the series (60). Doing that, the Eq. (58) 
becomes: 



/(*) = E 2 ~W^ I d ^ £ ^^ 2 l n s s ^- n ^^ /T cos (-nq,- • 0) . (61) 

j (particles) y A7T ) J n=l U 

Since the integrand function in Eq. (59) is periodical, the integration can be performed in the 
interval [— 7r,7r] instead of [0, 2tt]. The reason of this shift in the integration interval is that 
a considerable property of the function /(0) is the presence of a maximum at = 0, and, 
consequently, a very peaked maximum in the same point for the function exp |V/(0)] for large 
values of V. In this case, the saddle-point approximation can be used in calculating the integral 
(59). Therefore: 

i(particlcs) y Z7T ) J n=\ U 



2Jj + l r 3 7 ^e-v / ^^f/ T 



= /(0) - (q, • 0) 2 V ^^r / d 3 p — ^^ • (62) 

Let us define now a 3 x 3 real symmetric matrix A whose elements are: 

V(2J j + l) r A3 ^e-V^/r 

A kl = } — V^ — - / ^P [ ; Qi iQi k • (63) 

i(particles) ^ h ' (l±7 s 3 e ^^ i' ) l 

By using this definition, the Eq. (62) reads: 



Thus: 



/(0) ~ /(0) - • ^0 . (64) 

C(Q) - 7A3 exp [Vf(0)} J d 3 e iQ "* exp [-0 • A0] . (65) 

If V is large enough, the integration can be extended from [—it, n] to [— oo, oo] without affecting 
significantly the final result. Hence: 
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C(Q) " (2^ eXp[V/( ° )] Vd^A eXP[ "I QA " Q] • (66) 

Now we are able to write an approximate expression of chemical factors in the Eq. (22): 

Z(Q) = exp[-lQA-iQ] = GXp[ 2 QA ^^~^ ^ (6?) 

By using this approximation, the average multiplicity of primary hadrons (29) in the Boltzmann 
limit can now be written as: 



(„,)) = (2J 1 + l)-^—Y,> Ji 3 P e'^+^ /T e« A "".' 2 e-* 4 "'"/ 4 . (68) 



To summarize, in the large volume limit, chemical factors reduce to a product of two factors: 
the first corresponds to a traditional chemical potential whereas the second does not have a 
corresponding grand-canonical quantity; its presence is ultimately due to internal (i.e. quantum 
numbers) conservation laws in a finite system. Since: 

lim A" 1 = (69) 

the additional suppression factor exp[— qjA _1 qj/4] is negligible in the proper thermodynamic 
limit provided that vectors q^ are finite: the grand-canonical formalism is recovered. 

C Heavy flavoured hadrons production 

As shown in Sect. 4, the average multiplicity of primary charmed hadrons in events in which 
one cc pair is created owing to a hard QCD process must be calculated with the usual Eq. (28) 
in which the partition function Z is (see Eq. (33)): 

Z = Z 1 (Q°)-Z 2 (Q°,0). (70) 

The function Z\ can be written in the very same fashion as in Eq. (31): 

^(Q°) - 7T15 / *** eiQ0>ex PE^7?e-^>+^-^ fd 3 p loga-e-V^^- 1 ^)- 1 ] , 
(2nf J . j=1 (2vrf J 

(71) 
while the function Z-i can be worked out according to the same procedure depicted for the 
function Z 2 in Eqs. (26), (27) for the leading baryon effect: 

Z 2 (Q°,K) = -±- j d 5 e iC ^ j W e iK * 



x 



(27T) 

3 V 



exp E zrfe-***-™* + E7^/d 3 P log (1 - e -V^^/ T -^)^] , 

(72) 
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where the second sum in the exponentials in both Eqs. (71) and (72) runs over the three pion 
states and \Cj\ in Eq. (72) is the absolute value of j th hadron's charm. 

Henceforth, we denote by Q° and qj, q^, q^ three-dimensional vectors having as components 
electric charge, baryon number and strangeness, while charm and beauty will be explicitely 
written down. By using this notation, the average multiplicity of a charmed hadron with 
Cj = 1 turns out to be (the Boltzmann limit holds, cf. Eq. (29)): 

//„ vv _ 7 ^i(Q°-q J ,-l,0)-^(Q -q J , -1,0,-1) m) 

W J Z 1 (QO,0,0)-Z 2 (QO, 0,0,0) • [{6) 

Since the z functions of heavy flavoured hadrons are < 1, as shown in Sect. 4, a power 
expansion in the z/s of all charmed and anti-charmed hadrons can be performed from Zj = 
in the integrands of Eqs. (71) and (72), that is: 

expE Y/zje-^} ~ 1 + £ Y/zje-^ + W 7?7?**;e~ ,(,u+qi) '* ( 74 ) 

3 J i,3 

for Eq. (71) and 

exp [J2 Y s j Zje-^-W^} ~ 1 + £ zje-^- ; |c ^ + - £ T^'Zi^e -1 ^*^ -21 |c ^ (75) 

3 3 i,j 

for Eq. (72). Furthermore, the z functions of the bottomed hadrons can be neglected as they 
are <C 1 as well and beauty in Eq. (73) is always set to zero. 

Those expansions permit carrying out integrations in the variables ip, 04 and 05 in Eqs. (71) 
and (72). Thus: 

Zi(Q° - q„ -1, 0) ~ 5>?*iC(Q° - q, - Qi) , (76) 

i 

where the sum runs over the anti-charmed hadrons as the integration in 4 of terms associated 
to charmed hadrons yields zero. The ( function on the right-hand side is the same as in Eq. (32). 
Moreover: 

Zx(Q , 0, 0) ~ C(Q°) + E 7?7?***C(Q° - <b - q*) , (77) 

i.k 

where the index % runs over all charmed hadrons and index k over all anti-charmed hadrons. 
Owing to the presence of the absolute value of charm in the exponential exp[i |C J |'0] (\Cj\ = 1) 
in its integrand function, the function Z 2 (Q°, C, B, K) vanishes if K < and yields K th -order 
terms of the power expansion in Zj if K > (see Eq. (72)). Therefore: 

Z 2 (Q°,0,0,0) = C(Q°) (78) 

and 

Z 2 (Q°,-1,0,-1)=0. (79) 

Finally, inserting Eqs. (76), (77), (78) and (79) in Eq. (73) one gets: 

tin-)) = l^z- S<Tff*C(Q°-q J --tt) (80) 
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where the indices j, k label charmed hadrons and % labels anti-charmed hadrons. From previous 
equation it results that the overall number of primary charmed hadrons is 1, as it must be if 
c quark production from fragmentation is negligible. The average multiplicity of anti-charmed 
hadrons is of course equal to charmed hadrons one. The same formula (80) holds for the average 
multiplicity of bottomed hadrons in events with a perturbatively generated bb pair. If leading 
baryon effect is taken into account, the formula (80) gets more complicated, but the procedure 
is essentially the same. 

It is clear that a possible charm or beauty suppression parameter 7 C or 7&, introduced by 
analogy with strangeness suppression parameter 7 S , would not be revealed from the study of 
heavy flavoured hadron production because a single factor multiplying all Zj functions would 
cancel from the ratio in the right-hand side of Eq. (80). 

D On the definition of pressure 

In Sect. 4, we have dealt with pressure as a single well-defined quantity for the whole system 
of hadron gas fireballs. However, since the system has local collective flows, the definition of a 
single pressure is not a trivial one. We will now show that the best (as well as the most natural) 
definition is: 

p = T ^r> ( 81 ) 

where Z is the global partition function (see Eqs. (18)-(20)) and V the global volume defined 
in Sect. 2. 

If the temperatures and 7 S parameters of the fireballs are the same, as we have assumed through- 
out, the global partition function depends on the single fireball volumes only through the sum 
J2iL\ Vi, as shown in Eq. (20). Therefore we can replace the derivative in Eq. (81) with: 

' =T %v-m- (82) 

In order to develop this equation, we can use the expression (15) of the global partition function 
and write: 

T N f) N 

P=vE%>ogj; ^. Q -o )QO II % (9?) > (83) 

which is equal to: 

rp N o N 

P = ^E w ^ ■ ■ ■ > Q")7*7 lQ g II Z j(Q°j) ( 84 ) 

V i=l Q?,...,Q^ UVi 3=1 

by using the weights defined in Eq. (13). It is now possible to expand the derivative in Eq. (84): 

a N N a q 

^logll Zj(%) = Ysgy^zAQ*) = ^iog^(Q°) ; (85) 

i 3=1 3=1 % % 

the last equality is due to the dependence of Zj(Q®) only on the volume Vj. 
We can now write the pressure as: 
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N V 8 

P = TY,w(Q° 1 ,...,Q%)Y,V7w. logZ ^- ^ 

Q!,...,QS, i=1 l 

The expression TO log Zi(Q®)/dVi in the above equation is the pressure Pi(Q°) of the i 4h fireball, 
which is a well-defined one, as the fireball is a system at complete thermal and mechanical 
equilibrium by definition. Therefore the global pressure turns out to be: 

N y. 

p= mQ?,...,C&)£^(Q°). (87) 

The last equation now makes it clear that the definition (81) is the most natural definition of 
pressure for two reasons: 

1. for a given event, in which a specific configuration of fireball quantum numbers {Q?, • • • , Q%} 
is created, the global pressure is the average of the pressures of single fireballs, weighted 
by their extension through the factor Vi/V; 

2. in general, the global pressure is the average over all possible configurations of fireball 
quantum numbers according to their probabilities of occurrence if (Q?, . . . , Qjv). 

E Average multiplicities and volume fluctuations 

We want to calculate the leading correction to the formula (38) for particle average multiplicities 
in presence of volume fluctuations affecting chemical factors Z(Q° — nc\j)/Z(Q ). According 
to Eqs. (35), (36): 



„ oo 

« % » =/ dV P (V) £(Tl) n+1 7r^(n)£(n) , 

n=\ 

(88) 



where /,(„) = Z(Q° — nqj)/Z(Q°). If the volume fluctuations are not too large, one can expand 
the chemical factors fj< n \ around the mean volume V up to the first order term: 

f S (n)(y) ~ f Kn) (V) + f j{n) (V)(V - V) , (89) 

so that: 



((nj)) = Y,^ l ) n+1 lT 3 Vi J(n) f j(n) {V)+ 1 : s ^ j{n) f j{n) {V) / dVp(V)V(V-V) . (90) 

n=l J 

The first term in the series above gives rise to the formula (38), whilst the second term can be 
written as: 

i7^Vi 3{n) Vf ]{n) (7) = ^ir>z jin) (7)vf j{n) (7) , (91) 

o 

where D is the dispersion of the distribution p(V). This term is then about a factor D 2 /V 
smaller than the leading term. 
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Figure captions 

Figure 1 Behaviour of the global partition function Z as a function of electric charge, baryon 
number and strangeness, keeping all remaining quantum numbers set to zero, for T = 170 
MeV, V = 20 fm 3 and 7, = 0.5 

Figure 2 Behaviour of the non-strange baryon chemical factor Z(0, 1, 0, 0, 0)/Z(0, 0, 0, 0, 0) as a 
function of volume for different values of the temperature and a fixed value 7 S = 0.5 for 
the strangeness suppression parameter 

Figure 3 Results of hadron multiplicity fits for pp collisions at \/s = 19.4, 23.8 and 26 GeV. 
Experimental average multiplicities are plotted versus calculated ones. The dashed lines 
are the quadrant bisectors: well fitted points tend to lie on these lines 

Figure 4 Residual distributions of hadron multiplicity fits for pp collisions at y/s = 19.4, 23.8 and 
26 GeV 

Figure 5 Results of hadron multiplicity fit for pp collisions at y/s = 27 A GeV. Top: the exper- 
imental average multiplicities are plotted versus the calculated ones. The dashed line 
is the quadrant bisector; well fitted points tend to lie on this line. Bottom: residual 
distributions 

Figure 6 Results of hadron multiplicity fit for pp collisions at yfs = 200, 546 and 900 GeV. Ex- 
perimental average multiplicities are plotted versus calculated ones. The dashed lines are 
the quadrant bisectors: well fitted points tend to lie on these lines 

Figure 7 Residual distributions of hadron multiplicity fits for pp collisions at \/s = 200, 546 and 
900 GeV 

Figure 8 Freeze-out temperature values found by fitting hadron abundances in pp, pp and e + e~ collisions 



33| as a function of centre of mass energy; they are consistent with a constant value over 
an energy range of about two orders of magnitude. The error bars within horizontal ticks 
at \/s = 27 '.4 GeV pp collisions and at \fs = 91.2 GeV e + e~ collisions are the fit errors; 
the overall error bars are the sum in quadrature of the fit error and the systematic error 
related to data set variation (see text) 

Figure 9 Strangeness suppression parameters 7 S found by fitting hadron abundances in pp, pp and 



e + e~ collisions |33] as a function of centre of mass energy. A slow rise of 7 S from 19 to 900 
GeV in hadronic collisions may be inferred. At equal centre of mass energy, 7 S appears to 
be definitely lower in pp and pp collisions than in e + e~ collisions. The error bars within 
horizontal ticks at \fs = 27 A GeV pp collisions and at \/s = 91.2 GeV e + e~collisions are 
the fit errors; the overall error bars are the sum in quadrature of the fit error and the 
systematic error related to data set variation (see text) 

Figure 10 Local density and pressure of the hadron gas at freeze-out derived from the fitted param- 
eters as a function of centre of mass energy 

Figure 11 Dependence of fitted parameters and primary average multiplicities (top) on the mass 
cut-off in the hadron mass spectrum for pp collisions at \/s = 900 GeV 
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